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Abstract 

The contribution to the nuclear transverse response function Rt arising from two particle-two 
hole (2p-2h) states excited through the action of electromagnetic meson exchange currents (MEC) is 
computed in a fully relativistic framework. The MEC considered are those carried by the pion and 
by A degrees of freedom, the latter being viewed as a virtual nucleonic resonance. The calculation 
is performed in the relativistic Fermi gas model in which Lorentz covariance can be maintained. 
All 2p-2h many-body diagrams containing two pionic lines that contribute to Rt are taken into 
account and the relative impact of the various components of the MEC on R^ is addressed. The 
non-relativistic limit of the MEC contributions is also discussed and compared with the relativistic 
results to explore the role played by relativity in obtaining the 2p-2h nuclear response. 
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I. INTRODUCTION 



Two particle-two hole states in nuclei can be excited through the action of two-body 
electromagnetic (EM) currents that enter via Noether's theorem and the requirement of 
minimal coupling once an effective Lagrangian (embodying baryonic and mesonic degrees 
of freedom) is invoked. The currents reflect the dual role played by mesons in nuclear 
structure studies: mesons carry both the force which accounts for the binding of nuclei and 
the currents (in particular the EM current) which respond to external fields impinging on the 
nucleus. Accordingly, the natural classification into MEC and correlation currents follows. 
In particular, both substantially contribute to the inclusive, inelastic scattering of electrons 
from nuclei in the excitation energy domain extending through the quasielastic peak (QEP) 
to high inelasticity, the region we shall explore in the present paper. 

A long-standing issue of hadronic modeling is the difficulty of treating currents and 
interactions consistently, that is, of respecting gauge invariance by fulfilling the continuity 
equation. The model employed in the present study, the relativistic Fermi gas (RFG) model, 
has special advantages over most others: gauge invariance can be respected and fully rela- 
tivistic modeling, while not easy, can be undertaken. Given the focus of modern electron 
scattering experiments on kinematics where the energies and momenta transferred to the 
nucleus are large, and thus where relativity is expected to be important, one has strong 
motivation to explore such models in which fundamental symmetries can be maintained, 
even at the expense of accepting their obvious dynamical limitations. Given the necessity 
to break the problem down into tractable pieces, in the present work we forgo the issue of 
gauge invariance and focus on a fully relativistic description of the MEC contribution to the 
EM excitation of the 2p-2h states of the RFG. This is already a non-trivial computational 
problem and thus the issue of how gauge invariance can be maintained (which is presently 
also being explored in other work) is left for future presentation. 

Furthermore, we shall limit our attention in this work to the roles played by the pion 
and A in obtaining the MEC contributions specifically to the transverse response Rt- The 
relevance of the vr and A in for the physics of the quasielastic regime is well-established 
and, moreover, we wish to compare our results with the only existmg fully relativistic com- 
putation of the 2p-2h MEC contribution to Rt, namely that of 1] (quoted as DBT in the 
following). That this is to our knowledge the only previous study of this type is not surpris- 
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ing. Indeed, to carry out our project it has been necessary to compute a very large number 
of terms. Specifically, to get the direct contribution (see below) we have had to compute 
about 3000 terms, whereas to get the exchange contribution (which, as we shall see, comes 
into play only via the A-isobar) it has been necessary to compute more than 100,000 terms. 
The analytic manipulation of the traces of the relativistic currents has been done using the 
algebraic computer program FORM Q], and the numerical computation of the individual 
relativistic contributions requires up to seven- dimensional Monte Carlo integrations. 

To make direct comparison between our results and those of DBT we have used exactly 
the same form factors and A-width of the latter, although more modern versions of these 
quantities can straightforwardly be incorporated — in future work we will update our pre- 
dictions by doing so. In order to appreciate the importance of Lorentz covariance, we have 
calculated Rt for the RFG not only fully relativistically, but also in the non-relativistic limit. 



For the latter a few calculations are available, specifically those carried out in [l|, |3|,y,|5|. As 
we shall see, in the current work some important differences between the DBT results and our 
results in the high energy domain are found, notably in the non-relativistic approximations 
obtained in the two studies. 



II. TWO-BODY MESON-EXCHANGE CURRENTS IN FREE SPACE 

In Fig. n we display the free-space two-body isovector MEC entering in our calculation. 
Using the labelling in the figure and defining the four-momenta ki = p[ —pi and k2 = P2 —P2 
(the four-momentum carried by the virtual photon is then q = —ki — ^2) their relativistic 
expression is^ 

j;(A;i, k^) = -^l-llmfj^i^r^^) X r(2))3n(A;i)(i)n(A;2)(2)(fc2 - k,r (1) 
for the pion-in-fiight current (diagram (g)) and 

J^h, k,) = W^(^a) ^ ^(2))^ [n(fc2)(2)(7V)(i) - miMi'i'h)] (2) 

for the seagull (or contact) current (diagrams (e) and (f)). In the above 

n(A:)» = jr^, (3) 
k^ — fi^ 



^ In this section the current "operators" should actuaUy be regarded as the spin-isospin operator parts of 
the corresponding matrix elements between the two-particle states specified in the diagrams of Fig. ^ 
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FIG. 1: The two-body meson exchange currents in free-space. The thick hnes in the diagrams (a) 
to (d) represent the A propagation. 

fijt and M are the pion and nucleon masses, respectively, V is the (large) volume enclosing the 
Fermi gas and fnNN WInnI^'^ — 0.08) the pseudo- vector pion-nucleon coupling constant. 
The other coupling constants are /^^r^ = 1 and f-yj^NN = fnNN] for simplicity, here we omit 
the form factors (but they are included in the calculations of the response). The index [i) 
attached to the vertex operators distinguishes between the two interacting nucleons. 
The A current is derived from the Peccei pseudo- vector Lagrangian and reads 



J^(/ci, k2 



1 fnNNfnNAf' 



-yNA 



2M^^l 

(2) : h-rW 



— T. 



(2) 



^(t« 



X r(^))3 



J(l)(P<^'^2,g)75j 



(1) 



+ '^3 +3 



n(A;2)(2) + (l-2) .(4) 



It corresponds to diagrams (a)-(d) oi Fig. ^ In the above equations denotes the isobar 
mass, pa=Pi — Pb = Pi + Q and fnNA = 0.54, /^ata = 5 yield the strength of the coupling 
of the A to the EM and pionic fields, respectively. 

In Eq. (0} the following definitions have been introduced: 

1 



and 



J(a)/.(p, k, q) = {Akfs - ^-ff3)S^"'{p, Ma)- {-ijl^ + Qt,!^) 



j{b)^,{p, = ^ {-ipilt^ + ?m7/3) '^^^(p, Ma)(4A;^ - 7^^); 



(5a) 



(5b) 



the analogous expressions associated with the diagrams (c) and (d) are simply obtained 
through the interchange (1 ^ 2) as indicated in Eq. which also implies pa Pc and 

Pb ^ Pd {Pc = P2- q, Pd = p'2 + <?)• 

Furthermore, the Rarita-Schwinger (RS) A propagator, namely 

is used (we employ here the metric a^b^ = oPlP — ah). 

In the Appendix, the explicit expressions of Eqs. (jSt^-b) are reported. 



The non-relativistic limit for the currents 



The above given currents are the basic ingredients required in modeling Rnp- Since one 
of our goals is to assess the quality of a given non-relativistic approximation for i?^, we also 
need the corresponding matrix elements of the currents in that limit. Under the assumption 
that the momenta of the nucleons involved in the response are small in comparison to their 
rest mass (which may or may not be valid) one obtains^ 



J/(fci, ^2) 
Js{ki,k2) 



l4 



\/2 

(2M)2 f-nNNf-^TTNN 



kl + ^^l kl + ^^l 



X 



2 4ri(5^^ _t-2<5t{ -ri<5r^ 



^ Xl'^f^-k2Xa2 Xl'^(^-klXai ^ 

k^ + kf + fi^ "2 



and 

JA(fci,fc2) 



-A6r,-r2Sri,-TiSr^,-r2xi'^[Q X (^2 X O")] X<xi 



T2- 



(7) 



(8) 



(9) 



For brevity, the momentum labels implicitly include the spin (cr) and isospin (r) quantum numbers. 
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where 



, _ . ^ ^ 3Ma + 3Mi _ ^ 6Mi - 3MMa + 

(Mi-M2) 3Mi(MA-M) ^ 



and 

r2 



^ 2MA + 3M+|g^ 3Ml _, 6Mi + 3MMA-M2 

(Mi-M2) 3Mi(MA-M) ■ ^ ^ 



It should be noted, however, that this is not the onlypossible choice for a non-relativistic 



A current. For instance the one employed in P] and 



41 indeed has the same structure as 



Eq. Q, but occurs with coefficients given by the following expressions 

(2Ma + M) (2MA + 3M) 

^-^Mi-M2 Mi-M2 ' ^^^^ 

which are obtained by disregarding the last two terms in the numerator of Eq. (fTUk -b). Note 
that these terms have opposite sign and affect the magnitude of the coefficients A and B, 
in opposite direction, by less than 10%^. The non-relativistic A-current employed in DBT 
is again given by Eq. 0, however with coefficients 

A = , ^ ^ ^ ^ and B = 2A, (12) 

3(Ma - M) ' ^ ^ 

derived by setting M = Ma in the numerator of the right hand side of Eqs. (jlUb -b). The 
2p-2h response function is stable with respect to these different definitions of the coefficients 
A and B (the changes amounting to a few percent). 



Obviously, the expressions in Eqs. (|1HI can be obtained by using the simphfied form for the RS A 
propagator, S^''{p, Ma) = {^ + Ma)/ {p^ — M^)g^'^ , the other terms giving contributions which are either 
zero or proportional to 1/Ma and 1/M'^. Note that the second term, namely —■^^■^"'/S, never contributes 
to EM processes. 
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III. THE RESPONSE FUNCTION 



We turn now to the 2p-2h transverse response function, which can be obtained through 
the hadronic tensor. The latter is defined according to^ 

x{F\Ji\p[ApiP-2){p[p'2PW2MF), (13) 



where \F) represents the Fermi sphere, Ep = \/p'^ + M^, J = + J/^ = Jf + + J/^ and 
the factor 1/4 accounts for the identity of the 2 particles and 2 holes in the final states. 

Being concerned with Rt in the present work, actually we only need the spatial compo- 
nents of W^u (in fact only the fi = u = 1 and 2 components in the combination 11 + 22, 
since the unpolarized transverse response is the focus). That is, we have 

D f ^ (2vr)^\/ f dp[ dp'^ dpi dp2 , , 

E (^^^- - ^) [jKpi-PuP2-P2)JM-PuP2-P2) 

ar i,j=l T / 

-4(P'l - PUP2 - P2)Jj{p'l - P2,P2 - Pi) • (14) 

This expression explicitly displays the direct and exchange contribution stemming from the 
antisymmetrization of the 2p-2h final states in Eq. (fT^ . In Eq. (fT^ the Jj are the matrix 
elements of the currents in Eqs. (^, (0) and (0} in momentum space. Note that we define 
the excitation energy as u; = — go- 

By eliminating one of the integrations via the S^^^ , a simpler expression for the transverse 
response emerges: 

RT{q,u;) = J dp[dp^dp29i\p[\-kp)9{\p'2\-kp)eikF~\p,\)9{kp-\p2\) 

x6\u - (Ep, + Ep, - Ep, - Ep,)] V V ^ ( 5i, - 

^ ^ "^'^^ ^ tj^,^^Ep,Ep.Ep,EpA' Q'J 

jj{ku fca) fcs) - 4(^1, k2)J,{k[, k'2)] , (15) 



X 



^ In the previous sections, the nucleon momenta pi and p[ were arbitrary. From now on, we will indicate 
with Pi and p'^, respectively, momenta below (holes) and above (particles) the Fermi momentum kp. The 
pion momenta have, as before, ki = p[ — pi. 
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where now ip'^ = —q + Pi + P2 — Pi- The exchange term in Eq. ()15p depends on the pionic 
momenta k[ defined as k[ = p[ — P2, = p'2 — Pi- 

By exploiting the remaining 5-function and with some algebra, the integral defining the 
transverse response can be reduced to seven dimensions. In previous work where only the 
direct term and the non-relativistic limit were considered, it turned out to be possible to 
reduce the integration to a bi-dimensional one; but now for the relativistic expression, this 
is no longer possible, even for the direct term. 

As anticipated in the Introduction the traces of the relativistic currents typically generate 
a huge number of terms, so that it is not practical to report explicitly the expressions used 
in the numerical calculations of the relativistic responses. However, in the Appendix we 
give the formulae for the direct pion and pion/A contributions, the only ones having sizes 
suitable for publication. 

Numerical integrations have been performed using Monte Carlo techniques, varying the 
sample size until a standard deviation better than 1% is reached. The number of configu- 
rations required for such an accuracy is 10^ -r- 10^. A more stringent test of the accuracy in 
the numerical calculations can be obtained in the case of the direct non-relativistic contribu- 
tion, since this can be reduced to a two-dimensional integral [3] and done through standard 
quadrature. Also here the agreement is at the 1% level. 

The non-relativistic limit 

Before addressing the issue of the numerical evaluation of the fully relativistic problem, 
below we report the non-relativistic limit of the integrands in Eq. ()15|) . to be referred to 
as lZ^{ki,k2;q,uj) and 7l^{ki, k2; k[, k2; q,u!), for the direct and exchange contributions, 
respectively. This will allow us to connect the various pieces contributing to the transverse 
response to specific Feynman diagrams and to ascertain from whence the major contribu- 
tions arise. While the relativistic expressions for the transverse response function are quite 
cumbersome, the non-relativistic limit provides a relatively simple and controllable envi- 
ronment in which to check the correctness of the calculations. We also remark that our 
expression for will t urn out to differ somewhat from the one derived in DBT (no explicit 
formula is given there for the exchange part). 

For the purely pionic contribution and for the interference between the pionic and A 
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(a) (b) (c) (d) 

FIG. 2: The direct pionic contributions to the MEC 2p-2h response function. 




(a) (b) (c) (d) (e) if) 



FIG. 3: The direct pionic/A interference contributions to the MEC 2p-2h response function. 




FIG. 4: The direct A contributions to the MEC 2p-2h response function. 
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current we get, respectively, 



7^?^(fcl,fc2;g,^) 



(2M)' 



err ij 



mj 



(fc? + /i2)2(fc2+^2)2 



^2 ^2 

fnNNflTrNNf- 



U2 pi 



+ 



'yTVTV 



and 



k\rj.k\F^jsfj^{k\)FT^]^]^(k\) ^ k\j,k\F^p^j^{k'i,)FT^]^]^(k\] 



{kl + ^^lnkl + ^^l) 



{kl + ^^lm + ^^l 



2^2 



(16) 



{2My 



EE * 



err 2j 



X 



/^Af AT fnNAf-yNA 



p2 



X 



firNN f'ymr 



ttNN 



{kf)F^NN{H)F^NA{kt: 



2^ k^k-^^ 



{kl + ^^lnkl + ^^l] 



f-fnNN 



F^^AkDF.NAikl]^''''^^'-''^ 



(kj + fil) 



{kl)F^NNikl)F^NAikl) 



IT 



{ki+ixim+i^D 



+ (1^2) 



(17) 



From the different coupling constants in Eqs. (fTBj) and (fTTj) it is easy to identify the con- 
tributions arising from the pion-in-flight and seagull currents. The topologically different 
diagrams are shown in Figs. |21 and |21 
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Finally, for the contribution to TZ^ of the A current alone, one has 



err 

^fTrNNfivNAf-yNA jp2 ( 2\ 2 

9MX ""^^^^ 



2 /I 2\ 2 /I 2\ / 2 ^k-^rpj^-^ 2 H~ 

X < ( S + ^2)2 + ^ (fc2 + ^2)2 



+ (1 ^ 2) 1^, (18) 



where the first two terms on the right-hand side correspond to the diagrams (a)-(c) of Fig.|5 
and the last one to the diagrams (d)-(f). In this case six distinct diagrams contribute. 

In Eqs. f|TT)jl . (fTTjl and (fT^ ki and kx indicate the longitudinal and transverse components 
of the vector k with respect to the direction fixed by q. Furthermore, in the appropriate 
places, the hadronic monopole form factors 

A 2 2 

F^NN{k^) = (19a) 
F^NAik') = (19b) 

and the EM ones 

F^NN{q^) = (i_g2/A2^)2 ' (19c) 

F,NA{q') = F,Mq') (1 - ||) ' ~ S ' ^^^"^^ 

have been introduced. In the non-relativistic expressions the hadronic form factors have been 
taken in the static limit. The cut-offs have been chosen as in DBT, namely A^r = 1300 MeV, 
Kna = 1150 MeV, Al = 0.71 GeV^, A2 = M + Ma and A| = 3.5 GeV^. This choice clearly 
makes it possible a direct comparison between our results for and those of DBT. 

For completeness, we give also the formulae of the (smaller) exchange contributions to 
the integrand of Eq. (fT^. lZ^{ki, k^, k[, k'^] q,uj), in the non-relativistic limit. The purely 
pionic contribution is identically zero, as a consequence of charge conservation and of the 
fact that the photon does not couple to a neutral pion. For the interference between pion 
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FIG. 5: The exchange pionic/A interference contributions to the MEC 2p-2h response function. 
-<---« A---<---H «--->---« $.---<---f(. #---<---« «--->- 




FIG. 6: The exchange A contributions to the MEC 2p-2h response function. 



and A (Fig. E)) we have 



{2My 



Mm 



q 

Bq 



j^\k,, k2)jf{k[, k'2) + Jt\k,, k2)jnK, K 



ik2 X k'2)l 



{kl + ^,l){k'i+f^l) 



+ 



+ (1^2) 



8f^NNfj^NNfyNAf7rNA ^ \ " ^2)^? + " ^2)^2 " ' K){k2 " k'^) - {q ■ fc2)(fc2 
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2 



{q-k,)k'^ 



l + f^Dik'i + f^l 
{q ■ k',){k, ■ k',) , {q ■ k[)kl - {q ■ k2){k[ ■ k2) , 



{kl + f,l){k'i + f,l) {kf + f,l){kl + f,l 

The contribution of the A alone (Fig. E)) is instead 



(1^2) 



TL^ iki, k2] k[, k'^; q,uj) 



(2M)4 5Z5Z ' 



jp{k,M)JHKX. 



^f^NNfnNAf^NA ^2 ^ j^2 



+AB 



{ki ■ k[){kiT ■ k[^) _^ (fci ■ fc^)(fciT ■ fc2T) 



J2 



,/2 



+ (1^2) 



(20) 



+ (1^2) 



(21) 
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m (MeV) CO (MeV) 

FIG. 7: The non-relativistic transverse response function Rt^QjUj) according to the present cal- 
culation (solid) and to the one of DBT (dashed) for q = 550 MeV/c and q = 1140 MeV/c; in both 
instances kp = 1.3 fm^^ and €2 = 70 MeV. Only the direct contribution is shown. 

Equations (fTBj) . (fTTj) and (fTHj) could in principle be compared with Eq. (5.11) of DBT; 
however, the overall normalization of the latter is not correct, since its dimension is not 
consistent with its definition (namely of being the transverse part of the amplitude T given 
in Eq. (4.8) of DBT); moreover, the relative weights of the interference and A contributions 
with respect to the pionic one differ, in our calculations, by a factor 2 and 4, respectively, 
from those of Eq. (5.11) of DBT. These factors, however, are not able to explain the marked 
difference between our results and those in that paper. Note that although the authors 
of DBT write down exactly the same expressions as we do for the non-relativistic MEC 
currents, actually they state that the non-relativistic procedure to get their Eq. (5.11) is 
applied at the level of the hadronic tensor, that is by reducing the (cumbersome) exact 
relativistic response. 

In Fig. [7| we now compare our results with those of DBT, where the non-relativistic 
Rt (without the exchange contribution) is shown for q = 550 MeV/c (left) and for q = 
1140 MeV/c (right), with an atomic mass number of 56 and utilizing a Fermi momentum 
/ci? = 1.3 fm~^. The latter value is employed for the sake of comparison with DBT, although 
in fact it is more appropriate for heavier nuclei. 

It is clearly apparent in the figure that our predictions differ significantly from those of 
DBT: while the discrepancy is mild for moderate values of uj (roughly, those encompassing 
the QEP), it becomes striking at higher energies, namely in the region of the so-called dip 
and of the A-peak. Here our transverse response function in the proximity of the lightcone 
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turns out to be larger by about a factor two at g = 550 MeV/c and by over a factor three 
at g = 1140 MeV/c. 

Note that, in order to conform as closely as possible with the DBT approach, we have 
accounted for the initial state binding of the two holes by phenomenologically inserting 
a 70 MeV energy shift €2 in the energy conserving 5-function appearing in the response 
function^. In a non-relativistic framework this binding energy merely produces a shift of Rt 
toward higher u values. In the same spirit and to comply with gauge invariance, we have 
added (as in DBT) to the pion-in-flight current in Eq. (|7j) two terms, to be viewed as the 
coupling of the virtual photon to a fictitious particle of mass equal to the nNN cutoff, A^, 
in spite of the very minor role these terms play in shaping the transverse response. 

Furthermore, it is worth also remarking that our present calculations agree with the 
previous findings of 0, |^ . On the other hand, the computation of DBT does not agree with 
a previous one presented by the same authors in Q], notwithstanding that both of them 
employ the same set of values for the parameters and the kinematic conditions. Indeed, the 
former is lower than the latter by about 25%. 

In the light of these new results, we anticipate contributions from MEG effects in the 
region beyond the QEP of a magnitude not previously foreseen in most previous work. 
Because the studies of ^ I4 1 were applied to moderate momenta and to a restricted range of 
transferred energies (not exceeding 300 MeV, in fact) owing to their non-relativistic nature, 
it is only in the studies of DBT that suggestions for such large effects can be found. 

In the next section we shall explore whether or not the fully relativistic treatment and 
the inclusion of the exchange terms modify this finding. 



IV. RESULTS IN THE RELATIVISTIC REGIME 



Direct Contributions 

In Fig. ISlwe display our relativistic results both without and with the energy shift €2, for 
the same conditions as those of Fig. [7| The exchange contribution is neglected here as well as 
in all the other figures of this subsection and considered separately below. In Fig. IHlwe also 

^ This amounts to replacing Ep-^ + Ep^ by Ep-^ + Ep^ — €2 ■ 
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FIG. 8: The relativistic transverse response function RT{q,uj) at q = 550 MeV/c and q = 
1140 MeV/c calculated with e2 = 70 MeV (solid) and with £2 = (dot-dashed). Only the di- 
rect contribution is shown. The non-relativistic results are also displayed in order to shed light on 
the role of relativity in the response (dotted). For the sake of comparison the relativistic results 
obtained in DBT are displayed (dashed). In all instances kp = 1.3 fm^^. 

compare our relativistic predictions both with those of DBT and with our non-relativistic 
results. 

From the figure we observe the following. 

a) The fully relativistic calculation differs considerably under some circumstances from 

the fully non-relativistic approximation (i. e., the solid and dotted curves in Fig. |Sl 
respectively). Specifically relativity implies an increase of Rt over the non-relativistic 
results in the region up to the A peak. Specifically, the overall effects are small in the 
domain of the QEP, are modest in the dip region and are substantial in the region of 
the A-peak. In the region beyond the A peak, relativity instead yields a substantial 
reduction of the response with respect to the non-relativistic predictions. Of course, a 
hybrid approach can also be adopted in which spinor matrix elements and kinematics 
are kept non-relativistic, but the dynamic A propagator is used — this is discussed 
below. 

b) The present relativistic and DBT relativistic calculations are in essential agreement at 

q = 550 MeV/c and up to the A peak for q = 1140 MeV/c. 

c) In the dip following the A peak and above it up to the light cone, our transverse response 

is larger by about a factor two-to-three with respect to the DBT results. 
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m (MeV) CO (MeV) 

FIG. 9: Separate contributions to the transverse response function RriQ, in the non-relativistic 
limit at q = 550 MeV/c and q = 1140 MeV/c: pionic (dotted), pionic-A interference (dash-dotted), 
A (dashed) and total (solid); kp = 1.3 fm^^. The exchange contribution is disregarded here. 

d) At variance with the non-relativistic case, the binding energy not only shifts the rela- 

tivistic Rt toward higher energies, but also increases it. 

e) Roughly speaking two effects are especially important in determining the characteristic 

behavior seen in Fig. |S1 One is the peaking produced by the dynamic A propagation 
and the second is the basic phase space available to the two-nucleon ejection process. 
For instance, the strong rise of -Rt close to the light-cone at g = 1140 MeV/c comes 
from the growth of the available phase space. We have checked this by setting all of 
the current matrix elements entering in Eq. (fTH|) to unity and finding qualitatively the 
same behavior as in the full calculation. 

To gain deeper insight into the transverse response, we display its individual contributions 
in Fig. ini (for the non-relativistic case) and in Fig. El (for the relativistic one). In both 
instances the contribution of the A to Rt is seen to be overwhelming (although this A- 
dominance is mildly reduced by relativity), the more so at the larger momentum transfer. 
This feature strongly suggests the need for an appropriate treatment of the A degrees of 
freedom in obtaining EM responses over the whole kinematical regime under study. In 
particular the dynamical propagation of the latter appears to be essential, as discussed below. 
This in turn implies that a realistic accounting of the width of the A should be undertaken: 
here, again for ease of comparison with DBT, we have treated this problem as they did. 
Specifically, we have included the width through the replacement — ir{s)/2 in 

the energy denominator of Eq. (jH)), but not in the spinor factor. As far as T[s) is concerned, 
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FIG. 10: As in Fig. I^J but in the relativistic case. 

s being the square of the energy of the decaying A in its rest frame, we have adopted the 
same expression as in DBT. 

To assess the crucial role played in Rt by the A propagator more fully, in Fig. ^2 we 
display the relativistic transverse response (hence computed with the dynamical A propa- 
gator) together with those obtained by neglecting the frequency dependence of S^'^ {p, M^) 
(the width of the A being neglected as well) in the two possible choices of a static {po = 0) 
and of a constant A propagator. The extreme sensitivity of Rt to S^'^{p, M/\) is clearly 
shown by the results in the figure. 

One can understand even better the importance of treating correctly the A propagation 
by comparing the full relativistic results with a calculation employing the non-relativistic 
7A^A and ttA^A vertices, but retaining the relativistic, dynamical energy dependence in the 
denominator of the A propagator, p"^ — M"^ + iM^T . In other words, one can invoke a hybrid 
model in which everything is kept non-relativistic, but where the dynamic dependence in 
the A propagator is retained. The results are mixed. If one wishes to study only the region 
extending through the peak seen in Fig. ^Jout to the light-cone and wishes only to work at 
relatively small momentum transfers such as those in the figure, then this hybrid approach 
is reasonably good, incurring errors typically of 10 15%. 

However, this is not a general statement. First, if one addresses higher g- values, much 
larger relativistic effects are apparent and the hybrid approach is not very successful. For 
instance, at g = 2(3) GeV/c the typical error made at the peak is ~ 60(100)%. Secondly, 
if one focuses on the scaling region — as we shall in the follow-up study that is presently 
in progress — then even at relatively low values of q the effects of relativity (i. e., effects 
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FIG. 11: The relativistic transverse response function RT{q,uj) at q = 550 MeV/c and q = 
1140 MeV/c computed with various versions of the A propagator to ihustrate the sensitivity to 
the latter: exact propagator (sohd), static propagator (dashed) and constant propagator (dotted) 
(see text for the related definitions). In all instances €2 = 70 MeV and kp = 1.3 fm^^. 

other than the nature of the A propagator) are large. For instance, even at g = 1140 MeV/c 
one finds effects exceeding 100% in the scaling region. This is not apparent in the figures 
presented here, since the cross sections are so small in that region. However, this does not 
mean that the scaling region is irrelevant — quite the opposite, the total cross section is 
also very small (but measurable) and the 2p-2h MEC effects are far from negligible. The 
details will be presented in the near future in a separate paper. 

Next let us comment on several general features that characterize the response in the 
2p-2h sector in contradistinction with the Ip-lh case. In particular, the clean theoretical 
separation between the Ip-lh and the 2p-2h contributions holds only in the pure RFG frame- 
work. Here, moreover, the Ip-lh response is strictly confined to the well-known kinematical 
domain. In finite nuclei, or even in a correlated infinite Fermi system, neither of the above 
statements is valid any longer. We stick however to the pure RFG because this is the only 
model that allows us to maintain the fundamental principles of Lorentz covariance, gauge 
invariance and translational invariance. 

The 2p-2h response function resides on the entire spacelike domain of the [uj, g)-plane, 
in contrast to the Ip-lh case, which, as mentioned above, is constrained to a limited region 
- fo. .ece. f„..y rela«v.t. .ud.es of *e Ip-lh .espouse see fl fl 3. S.ce .o.e.u. 
conservation involves only the total momentum there is no limitation on the momentum of 
a particle taking part in the 2p-2h excitations and hence the energy of the latter will vary 
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over a broad range extending from lo = to the light-cone. 

In summary, naively the 2p-2h response function RT{q,uj) might not be expected to 
display any particular structure, but rather to appear as a broad background. As a matter 
of fact, as we have seen, RT{q,uj) indeed displays structure, but related not to the nuclear 
spectrum involved but to the role played by the A resonance. Furthermore, we have found 
that it grows with u, reflecting the associated growth of the phase space. Indeed the density 
per unit energy of the 2p-2h excitations goes as u^. 

Exchange Contributions 

n 

Considering the purely pionic MEC currents, it is well known 10] that the Ip-lh excited 
states of the RFG are reached only through the exchange term of the related MEC diagrams. 
In contrast the 2p-2h states are reached only through the direct term. The origin of these 
kinds of "selection rules" stems, in the first case, from the spin-isospin saturation of the 
RFG, in the second from the impossibility of fulfilling charge conservation in the exchange 
contribution for the diagrams displayed in Fig. |21 

The A-current, however, eludes the above constraints. It provides a direct contribution 
in the Ip-lh sector of the RFG spectrum, as well as an exchange contribution in the 2p-2h 
sector, both alone and through the interference with the pionic currents. 

The exchange contributions are included together with the (larger) direct ones in Fig. El 
Two features are immediately apparent from the figure. First, the impact of relativity on 
the exchange terms is as substantial as it is for the direct ones. Moreover, as for the direct 
case, the exchange response is first enhanced in the region of the A and then dampened by 
the relativistic effects as one approaches the light cone. 

Second, the relative magnitude of the exchange terms is seen to increase both with the 
energy transfer and with the momentum transfer. The fraction of the total relativistic 
response due to exchange goes from roughly 4% at low uj to 29% at the light cone for 
q = 550 MeV/c and from 7% to 39% for q = 1140 MeV/c. Thus, the exchange contribution, 
while not the dominant one in the kinematical situations we have explored, appears to us 
not only far from negligible, but also potentially interesting in connection with scaling (see 
below). 

Finally, as for the direct case, note that the present calculation of the non-relativistic 
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FIG. 12: The transverse response function Rxiq, uj) at q = 550 MeV/c and q = 1140 MeV/c includ- 
ing the exchange contributions: non-relativistic direct (positive dotted), non-relativistic exchange 
(negative dotted), non-relativistic total (light solid), relativistic direct (positive dashed), relativis- 
tic exchange (negative dashed) and relativistic total (heavy solid). In all instances €2 = 70 MeV 
and kp = 1.3 fm~^. 

exchange contribution agrees with the calculation at lower momenta reported in P| , whereas 
does not agree with that of DBT. The relativistic response, on the other hand, roughly agrees 
with DBT. 



V. CONCLUSIONS 



In this paper we have carried out a fully relativistic calculation of the MEC contribution 
of 2p-2h excitations to the transverse inclusive response within the context of the RFC 
These excitations are reached through the action of the MEC including the part associated 
with the A degrees of freedom. We have compared our results with those of DBT, the only 
existing relativistic study before the present one and we have pointed out similarities and 
differences between the two studies. To gauge the difficulty of this project it might help to 
stress again, as anticipated in the Introduction, that the number of terms we have computed 
to get Rt amounts to over 100, 000. 

Several new aspects contained in the present study bear highlighting. In particular, while 
the fully relativistic results obtained here roughly agree with those reported by DBT, some 
substantial differences occur, especially as the lightcone is approached. Furthermore, our 
non-relativistic results are completely in accord with our own previous work; however, they 
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are not in agreement with DBT in many circumstances. Indeed, the strongest disagreements 
occur at high inelasticity and suggest to us that one should exercise extreme caution when 
invoking traditional non-relativistic approximations in such a regime. 

Going beyond the extreme non-relativistic approximation [ll|, one can invoke various 



classes of hybrid mode 



model are retained 



ing where some, but not all, of the features of the fully relativistic 
isi [l4 1 . We have pursued such an approach in which the dynamic A 



propagator is kept, but otherwise everything is taken to be non-relativistic. The results are 
interesting: at modest q (say 1 GeV/c or less) and at high uj the hybrid approach is quite 
good, incurring only 10 -^ 15% error. However, at higher g-values, even in the high uj region, 
and at low cj, even for modest g, the hybrid approach is not very successful and relativistic 
effects are important. 

The high uj region, including the peak produced by the dynamic A propagator and the 
rise seen as the light-cone is approached when q is large, might pose some problems for 
ligh-energv photoreactions. The latter are related in the familiar way (see, for example, 
15l |) to the inclusive electron scattering transverse response by 

rv 

= 27r2-i?T(g,g), (22) 

UJ 



a being the fine-structure constant. 

In context, it is important to note that, while one cannot go 



rom these results directly 
Q 3), since the RFG 



to studies of (e,e'NN) or (7,NN) reactions (see, e. g., Refs. Il6 
is not well-suited to modeling such semi-inclusive processes, the understanding gained in 
the present work concerning relativistic effects, A propagator effects, etc., likely can have 
an impact, for instance via improved descriptions of current operators for use in other 
approaches. The A propagator provides just one example: it produces a prominent peak in 
the 2p-2h transverse response and apparently is an important contribution in the total cross 
section. It is also very dependent on the specifics assumed, such as what model is taken 
for the propagator, how medium effects are or are not incorporated, how ofF-shellness is 
handled, etc. In particular, in view of continuing discussions about the off-shell behavior of 
the A inside the nuclear medium and of our findings in this work, the 2p-2h transverse 
response may provide an ideal place in which to explore such issues. 

Clearly, this is not the intent of the present work. Here we wish only to make contact 
with the work of DBT, using exactly their ingredients and only in future work will we return 
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to address some of these interesting issues. 

In particular, in the near future our intent is to address the difficult issue of the gauge 
invariance as was explored in the Ip-lh sector in . This requires the introduction 

of 2p-2h correlation currents in addition to the MEC. Only when this task is accomplished 
will a consistent (namely, gauge invariant), Lorentz covariant and translational invariant 
description of the inclusive responses (both longitudinal and transverse) within the RFG 
framework be achieved. 

Furthermore, we are now in a position to study the degree to which these responses 
fulfill or violate y-scaling of first typ e (namely in the momentum transfer q) and of second 
type (namely in kp) — see jsof 0, 122]. An exploration of these issues is currently being 
undertaken for the 2p-2h MEC and will be presented in a forthcoming paper, which will also 
focus on the important theme of the interplay between the scaling violations arising from 
the MEC action in both Ip-lh and 2p-2h sectors. 
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APPENDIX 

For completeness, in this appendix we report some of the detailed formulae used for the 
relativistic calculation. 

We start by giving the explicit expressions for Eqs. (jSt^-b) which enter in the relativistic 



22 



A-current: 

-V^h^ii> - 4 (A; ■ p)p^ \^ q^p^ 

+3^ (4 {k ■ p) ri^ - p' h'i + 2 (g • p) h'^ 

-2p^ Ui) -A{k-p)q''^ + q^'p^ l/c^ 
+ (^7"^ - 2 -Sq^'^^ + Aq^ik-p) ij 

J(bMk,q) = 7o(i(l)(p,/i;,g))Vo (A. lb) 

Useful relations follow from momentum conservation: q^ + k'^ + k2 — 0, kiT + fe2r = and 
q — qL — —kiL — k2L- We also reiterate the definitions: pa=pi — q, Pb = Pi + q, Pc = P2 — q, 
Pd = P2 + ki=p[- Pi. 

Next, in the following formulae, we quote explicitly the integrand entering in the direct 
contributions to the transverse response, limiting ourselves to the full pionic contribution 
and to the pionic-A interference, since the pure A contribution corresponds to a far too 
lengthy expression to be reported here. We use the notation £ = Ep^Ep'^Ep^Ep'^. The 
hadronic form factors are not explicitly displayed. 

Pionic contribution: 

n^^ki, k2;q)^ (^^^ - ^) Wiki, k2) 

= i^E E J::^Hk„k2mk„k2) 

ar m=x,y 



^2 . .4 c { fnNNf' 

i 

'^fnNNf-ywNN 



{kl-^xl){kl-^^l) V M^{kl-^^l 



\2 



f'KNNf'y^NNf'yK-K 



(ki - fiinki - fii: 



+ (1^2) 



(A.2) 
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A-pion interference contribution: 



UT m=x,y 



+ fnNNf-rirNN 



fnNNfjTrTT 



III, ^ ' -' 



+ (1^2) 



k^oiS. 



AM{kl - lil) 



+(1^2) 



;(A.3) 



pure A contribution: 



ar m—x,y 
p2 f2 f2 



/2 £2 f2 
wNnIttNAJ^NA 



E 



m=x,y 



^ripup'i'^i)^riP2,P2,q) 
m-i4){ki-i^i) 



"^2 



(u2 2 ^2 { - ^).^(™) (^'^ ^2, q) {^1 + M)j^,) (p„ k2,q)} 
+Tr {i^, + M)j(-)(p,, ^2, qW - M)j}^){p,, k2, q)} 



-Tr {(A + M)i^)(p,, A;^, - M)j(-)(-p,, /c^, g)} +(1^2) 



^Tjpi^Pi, g)-^r(P2,P2; g) + ^ir{pi,Pi, q)^riP2,P2, g) 

3(A:f-/x^)(A:i-/z2) 



(A.4) 
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In the previous equations the following definitions have been introduced: 



KiP,P';Q) = Tr{{^-M){^' + M)j};^{p-q,-q-p'+p,q)} 
= Tr {{^ - M){f + M)3^^^{p^, q)] 
1 I 4g2 



pI-MI] 3Mi 

4M2 



p'"(p-g)-p"(p'-g) + 2p"^(p-p') 



pI 



4:{K-Pa)\ 



3Mi 



{^K ■ pa) [p'^ip' ■ q) - P"^{P ■ q) + ^p'^q^] - qV^ip^ + p'^) } 
^{K ■ Pa) {q ■ Pa) - pW + 2(? • Pa)' 

q^iK-PA) + {K ■ q){q-A) + 2{K ■ q){K-PA) 



+ 3Ma I " 

+4 {g^ [p'^{q ■ p) - p'^iq ■ p')] + 2(p • p') [p'^q^ + 2p'^{q ■ p') - 2p""(g • pa)] } 

+4M2 |4p'"^ [(g • pa) - P'^ip' ■ q)] - q'^ip'"' +p'^)} + SMMaC?" | , (A.5a) 
Tr{{^ + M) (/ - M) ip' + q,-q-p'+ p, q) } 

Tr { if -M){^ + M)j^^^ {p' + q^-q-p' + p^q)}^J^ (p', p.-q). (A.5b) 
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Sa{p,p';q) ^ Tr{{^-Mh^{f + M)j^^^{p-q,-q-p'+p,q)} 

m=x,y 

= Y.Tr{{i>- M)7"^(/ + M)j(-) (pa, q) ] 

m=x,y 

= I {4(A;^ ■ Pa) 2(g ■ pa) ({p'-a ) - (g ■ pl)) - ■ pl) 

-3(g -pOl^vr - Pa) - {q ■ Pa){K ■ P') + fcyrT " PT + " P')) 

-Pt-p'Aq -K) I + ^^|8(g ■PA)(/i;^ -Pa) 
+^j^\{q-p){K-p) - {q-p){K-p') - {q-K){p-p') -4:{pa-K){p-p') 



q^{K-p) - {q ■ K){p' -Pa) 



3Ma 
+Pt ■ Pt [pI - {p 



16M^ 
3Ma 



(q-K) ip-p')+PA -iq-p)iK-p') 



+2q\K ■ p') + 3(g • p)iK ■ p) - 4(g ■ p'){K ■ q) + 2(g • pa){K ■ Pa) 
-Pt + (q ■ P)) + Pt ■ p'Aq ■ Pa)] + 8M\^q^{K ■ p) ~ {q ■ K){p' ■ Pa) 

+iq ■ p'){K ■ Pa) - {q ■ pa){K ■ p') + Kt ■ Priq ■ p') - '^Kt ■ p'Aq ■ pa) 



+Pt 



q^ + 2{q-p') 



PtPt 



{q-p) + {q- p') 



iq-K 



+16Ma|(p ■ p') y + 2(g ■ p')] +Pt-p't + " p')] - PriQ " P)} 
+16MAM2|g2_2(g.p)| 



(A.6a) 



and 



Sb{p,p';q) ^ Yl Tr{{^ + M)j"^{f-M)jj;^{p' + q,-q-p' + p,q)}=Sa{p',p;-q). 

(A.6b) 



m=x,y 



In Eqs. ()A.5a|) and ()A.6a|) we used pa = p—q and /e^ = —q—p'+p to simplify the expressions. 
Eqs. frOb|) and (fA^ are consequences of the identities: j^^h){p,k,q) = j}^^i,){p, —k, —q) 
and Eq. (lAlbl) . 

The traces indicated in Eq. ()A.4|) involve thousands of terms and are not included here. 
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From Eqs. ()A.2|) - ()A.4|) the non-relativistic limit of Eq. p7p can be obtained by assuming 
that all the nucleonic momenta are much smaller than the nucleon rest mass. 
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